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Abstract. The isolated neutron star RX J0806.4-4123 was observed with XMM-Newton in November 2000. The 
data from the three EPIC instruments allowed us (i) to derive an improved X-ray position to an accuracy of 
2-3" , (ii) to accumulate the first medium-resolution soft X-ray spectra of high statistical quality and (iii) to find a 
candidate for the neutron star rotation period. Although this period of 11.3714 s is formally detected at a 3.5a level 
in the EPIC-pn data, the similar pulse profiles deduced from all three EPIC instruments increase the confidence 
that the period is real. The pulsed fraction of ~6% would then be the weakest X-ray flux modulation detected 
from dim isolated neutron stars. We fitted the X-ray spectra with blackbody and neutron star atmosphere models 
and discuss the results with respect to the brightness limit placed by optical images. The reduced size of the error 
circle on the X-ray position should allow deeper searches for an optical counterpart. 
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1. Introduction 

The ROSAT all-sky survey source RX J0806. 4-4123 
(1RXS J080623.0-412233) was sugge sted as candidate for 
an isolated neutron star (INS) by Haberl ct al. (1998) 



based on X-ray properties of this source similar to the es- 
tablished cases R XJ1856. 5-3754 (|Walter et al. 1996| ) and 
RXJ0720.4-3125 (|Haberl et al. 1997fl. While the latter 



two (the brightest in the sample of INSs discovered with 
ROSAT) were optically identified with my ~ 25.6 (Walter 
fc Matthews 1997|) and m B ~ 26 flMotch & Haberl 1998|; 



Kulkarni & van Kerkwijk 199?) blue objects, respectively, 
optical observations of RX J0806. 4-4123 were only able to 



emission has been detected yet from these objects. While 
little doubt is left that these objects are indeed INSs, the 
origin of their X-rays is still a matter of debate. If we 
see thermal emission from young cooling neutron stars, 
the observed blackbody temperatures in the range of 50- 
120 eV imply ages younger than ~ lfj 6 years for standard 
cooling curves. Older neutron stars may be re- heated by 
accretion of interstellar matter, producing X-ray luminosi- 
ties sufficiently high to be observable. However, the accre- 
tion rate strongly depends on the velocity of the neutron 
star relative to the interstellar medium, and the neutron 
star must rotate slowly enough to allow the matter to 



derive a limiting D magnitude mn > 24 for a possible opti- 



enter the magnetosphere (see the review of Treves et al 



cal cou nterpart ( Haberl et al. 1998 ) . High resolution X-ray 



spectr a obtained from RXJ1856. 5-3754 and RXJ0720.4- 



2000). The recent estimates of a pulse period derivative 
of (0.7 - 2.0) • 10" 11 



3125 h ave revealed no evidence for photospherie absorp- 

tion or emission features and are well modeled as black- 



s s from RBS 1223 (Hambaryan 



body radiation with little interstellar absorption (Paerels 



et al. 2001; Burwitz et al. 2001) 



Today seven ROSAT discovere d INSs (hereafter dim 



INS) are kno wn (see references in Motch 199S ; Zampieri 



et al. 200l|) with RX J0 720.4-3125 (rotation period 



P=8.39 s; Paberl et al. 1997|) , RX J0420.0-502 2 (P=22.7 s; 
Haberl et al. 1999| ) and RBS 1223 (P=5.16 s; [Hambaryan 
et al. 2002 ) showing pulsations in the X-ray flux. No radio 



ct al. 2002|) and - 10~ 14 s s" 1 from RX J0720. 4-3125 
(Zane et al. 2002) with inferred magnetic field strengths of 
~2- 10 14 G and ~ 10 13 G, respectively, favour the models 
invoking a young (or middle-aged) cooling neutron star 
rather than the accretion heating, at least for these pul- 
sating sources^. By virtue of this, the dim nearby neutron 
stars are suggested to be different from anomolous X-ray 
pulsars which have been prop osed to be powered either by 
superstrong magnetic fields (Thompson fc Duncan 1996 ) 
or accretion ( Mereghetti fc Stella 199q )^ 



Send offprint requests to: F. Haberl, e-mail: fwh@mpe.mpg.de 

* XMM-Newton is an ESA Science Mission with instruments 
and contributions directly funded by ESA Member states and 
the USA (NASA) 



1 [Kaplan et al. (20*02 ) derived a less stringent upper limit for 
P for RXJ0720.4-3125 by not including XMM-Newton data, 
but also their results support the interpretation that this source 
is a cooling neutron star. 
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Table 1. XMM-Newton EPIC observations of 
RXJ0806. 4-4123 on November 8, 2000. 



Camera 


Read-out 


Filter 


Observation 


Exp. 




Mode 




Start End (UT) 


[ks] 


MOS1/2 


FF, 2.6 s 


Thin 


13:43 18:47 


18.0 


pn 


FF, 73 ms 


Thin 


14:24 18:53 


15.6 



and are well modeled as blackbody radiation with little 
interstellar absorpt Following the spin history of the three 
known pulsars among the seven dim INSs and searching 
for pulsations in the others with new sensitive instruments 
is important for the understanding of the origin of their 
X-ray emission. Here we report on an XMM-Newton ob- 
servation of the poorly known INS RX J0806. 4-4123. We 
present an improved X-ray position, the X-ray spectrum 
and the detection of a candidate for the rotation period 
of the star. 



2. The XMM-Newton observation 



XMM-Newton fljanscn et al. 200l|) observed 



RXJ0806.4-4123 during the satellite revolution #168 on 
November 8, 2000. Here we report on the analysis of 
the data collected with the European Photon Imaging 
Camera s (EPICs) based o n MOS (EPIC-MO S1 and 



MOS2, Turner et al. 2001) and pn (EPIC-pn, Btruder 



et al. 2001) CCD detectors which are mounted behind 
the three X-ray telescopes ( Aschenbach et al. 2000| ). 
All cameras were operated in the Full-Frame mode 
(FF) providing data over a field of view of ~13' radius. 
Further details of the EPIC observations are summarized 
m Table |. The data were processed using the XMM- 
Newton analysis package SAS version 5.2 to produce the 
photon event files and tools from version 5.3 for further 
spatial, spectral and timing analysis. The event times 
were corrected to solar barycenter using the SAS task 
"barycen 1.13.1" (in which a timing problem present in 
older versions has been solved). 

2.1. The X-ray position of RXJ0806. 4-4123 

Haberl et al. (1998) derived an error radius of 11" 
(90% confidence) from the ROSAT PSPC observations of 
RX J0806. 4-4123. The error is dominated by systematic 
uncertainties in the attitude reconstruction which could 
not be reduced due to the lack of other X-ray sources with 
optical counterparts detected in the field of view to define 
a reference coordinate frame. In order to find field stars 
in the EPIC images, a standard source detection analy- 
sis using combined box sliding and maximum likelihood 
techniques as available in SAS, was performed simulta- 
neously in the three energy bands 0.3-1.0, 1.0-2.0 and 
2.0-7.5 keV. This yielded 50, 29 and 33 detections with 
existence likelihood larger than 10 in the pn, MOS1 and 
MOS2 images, respectively. Fig. |l| shows an EPIC-pn soft 



X-ray image. To identify possible optical counterparts to 
the X-ray sources, finding charts were produced using the 
DSS2 (red) image. To secure the identification, the fol- 
lowing criteria were applied: (i) only X-ray sources which 
were detected by at least two EPIC instruments were used, 
(ii) a stellar-like object was found in the X-ray error cir- 
cle, and (iii) both X-ray (spectral) and optical properties 
(R and B magnitudes as derived from the USNO A2.0 
catalogue produced by the US Naval Observatory) must 
be compatible with a star. This resulted in the identifica- 
tion of five X-ray sources with field stars as summarized in 
Table |[ Hardness ratios are defined as HR1 = (B+A)/(B- 
A) and HR2 = (C+B)/(C-B) with A, B and C as count 
rates in the energy bands 0.3-1.0, 1.0-2.0 and 2.0-7.5 keV, 
respectively. Coronal emission from stars is characterized 
by relatively soft X-ray spectra with the bulk of emis- 
sion below 2 keV. This results in negative values for both 
hardness ratios with HR1 depending on the amount of in- 
terstellar absorption (see also Haberl 2002 ). HR2 values 
of weak sources are often undetermined due to the lack 
of statistics in the 2.0-7.5 keV band. For comparison the 
X-ray properties of the target are also shown in Table ||. 

Using the SAS task "eposcorr" the X-ray source po- 
sitions were corrected to the USNO A2.0 catalogue ref- 
erence frame by minimizing the differences in X-ray to 
optical position for the five field stars. The uncertainties 
on the required coordinate shifts can be regarded as the 
remaining systematic error of the X-ray bore-sight and 
are 2.1", 3.1" and 2.5" for pn, MOS1 and MOS2, respec- 
tively. Fig. H shows the B-band image of the area a round 
RXJ0806. 4-4123 as obtained by |Haberl et al. (1998| ) with 
the three error circles derived from the EPIC instruments 
after the bore-sight correction. The corresponding posi- 
tion deduced from the pn image is RA = 08 h 06 m 23!47 
and Dec = -41°22'32"3 (J2000.0). The error circles ob- 
tained from the EPIC data agree with each other and are 
completely inside the ROSAT circle. This confirms that 
none of the surrounding stars is responsible for the X-ray 
emission from RX J0806.4-4123. 



2.2. The temporal behaviour of RX J0806.4-4123 

EPIC light curves of RX J0806. 4-4123 were extracted in 
the 0.12-1.2 keV energy band (using events with pixel 
patterns 0-12) which covers most of the X-ray spectrum 
(see the next Section). The light curves show no signifi- 
cant variations over the duration of the observation with 
average net count rates of 1.81 cts s _1 , 0.369 cts s _1 and 
0.405 cts s" 1 in the pn, MOS1 and MOS2 detectors, re- 
spectively (Fig. ||). To search for periodicities on short 
time scales, we used 29,549 events extracted from a circle 
with radius of 45" centered on the position of the target 
in the pn data. Of these selected events, about 97.5% were 
estimated to belong to the sour ce. As a periodicity te st we 
used the Rayleigh Z% method ( Buccheri et al. 1983 ) with 
one harmonic involved (n = 1) which is expected to be op- 
timal to search for rather smooth pulsations. The 0.073 s 
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Table 2. Field stars near RXJ0806. 4-4123 observed with EPIC pn. 



ML 1 


count rate 2 


HR1 


HR2 


USNO counterpart 


run 


771B 




[lO" 3 s" 1 ] 








[mag] 


[mag] 


518 


30.0±2.5 


-0.37±0.08 


-0.68±0.15 


U0450.05909836 


13.0 


14.0 


124 


8.5±1.1 


-0.22±0.13 




U0450_05911062 


15.0 


16.0 


74 


8.4±1.4 


-0.41±0.15 


-0.11±0.32 


U0450.05910649 


15.3 


16.5 


46 


3.7±0.8 


-0.61±0.19 




U0450.05889986 


17.3 


19.2 


28 


4.8±1.2 


-0.17±0.26 




U0450_05866926 


13.3 


14.5 


10 5 


1060±10 


-0.96±0.01 


-0.92±0.06 









1 Likelihood for existence with probability pml = 1 — exp(— ML) 

2 0.3-2.0 keV 



05910649 
m 

. .■ ™ , 



i 

05866926 



• - - # c > $m 

1 - a 4 | m.\ 1 4 . I " J- J ■ 




05911.062. y 
■ m. - 



05909836 



10' 



Fig. 1. EPIC-pn image in the 0.2-2.0 keV band. The 
image is smoothed with an adaptive intensity filter and 
out-of-time events (caused by photons detected during 
CCD-frame read-out) are subtracted on a statistical basis. 
Marked are the identified USNO A2.0 objects (see Table |) 
used to define the coordinate reference frame. 



time resolution of the pn data limits the periodicity search 
to frequencies below ~ 3 Hz. We performed the Z\ test 
in the A/ = 0.001 — 3 Hz frequency range with a step of 
5 • 10 -6 Hz. The oversampling by a factor of 10 (compared 
to the expected widths of 1/T sp of the Zf peaks, where 
T sp = 15.7 ks is the observational time span) was chosen 
to resolve separate Z\ peaks and to assure that no peak 
corresponding to a periodic signal is missing. As the vari- 
able Z\ has a probability density function equal to that of 




Fig. 2. B-band CCD image around the X-ray position of 
RXJ0806.4-4123 (from |Haberl et al. 1998j ). The circles 
represent the position uncertainties derived from pn (ra- 
dius 2.1"), MOS1 (3.1") and MOS2 (2.6") which include 
statistical and (dominating) systematic X-ray errors and 
a 0.25" error from the USNO A2.0 positions. 



Z\ > 2 log(JV) = 21.5 (peaks with Z\ < 21.5 are regarded 
as noise with a probability of 100%). This test resulted in 
only one peak higher than that lower limit, Z\ = 33.9 at 
the frequency /* = 0.087936 Hz (see Fig. |). The proba- 
bility to obtain by chance a peak of Z\ = 33.9 in N inde- 
pendent trials is p = N exp(-Zf/2) = 2.0 • 10~ 3 , which 
corresponds to a detection of pulsations at a confidence 
level of C = (1 -p) ■ 100% = 99.8%, or 3.1cr. Involving 
higher harmonics gives values Z\ = 42.9, Zf = 43.0 and 
43.2 at the same frequency /*. These numbers indi- 



zl 



cate that, according to the if-test suggested by de Jager 



X wit i 2 • n degrees of freedom (e. g., Bendat & Piersol et al. (1989), there is a statistically significant contribu- 



1986(^7 the probability to obtain a noise peak of a given 
Z\ height in one trial is exp(— Zf/2). The number of in- 
dependent trials in the chosen frequency range estimated 
as N — Af T sp = 47. 1 • 10 3 puts a lower limit on a Z\ value 
to be yielded by a periodic signal with a probability p > 0, 



tion from the second harmonic. A probability to obtain a 
noise peak of Z\ = 42.9 in one trial is only p = 1.1 ■ 10~ 8 
and in N independent trials is p = 5.1 • 10~ 4 , that in- 
creases the significance of the period detection up to a 
3.5cr level. We note this Z| peak is highest in the frequency 
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range chosen for the periodicity search (all other Z\ peaks 
are below 27.0 and produced by noise with probabilities 
p > 94%). We also performed the Z^-test adding 14,210 
source photons detected with two MOS cameras to the pn 
data. The increase of the total statistics by about 50% re- 
sults only in small increments in Z 2 values (Zf = 37.5 and 
Z\ = 43.9 at the same frequency /*) and significance of 
the detected signal (to a 3.6a level). That can be explaind 
by the low time resolution of the MOS data (see Table § 
which should lead to a strong smearing of the weak signal. 

To determine the pulsation frequency and its un- 
certainties more accurately, we applied the Odds-ratio 
method by Gregory fc Loredo (1996 ) based on the 



Bayesian formalism. The frequency-dependent Odds-ratio 
0(f) specifies how the data favour a periodic model of 
a frequency / over the unpulsed model. This ratio al- 
lows one to obtain the probability distribution function 
p(f) oc 0(f) f~ l for a signal to be pe riodic in a chosen 




frequency range (see Zavlin et al. 2000, for more details) 



Fig. U shows the dependence of O(f) in a range around 
the frequency /» resulting from the Z\ test. The maxi- 
mum value O max = 1.3 • 10 5 is at / = / = 0.087940 Hz 
(f -f* = 4-10~ 7 Hz < 1/T sp ). The uncertainties of f at 
the 68% and 90% confidence levels are 5f = 2.3 • 10~ 6 and 
10.5 • 10~ 6 Hz, respectively. The corresponding rotation 
period is Pq = 11.37139 s with uncertainties of 3.0 • 10 -4 
and 13.6 • 10~ 4 s at the 68% and 90% levels, respectively. 

The folded light curve extracted at / = /o from the 
pn data (see Fig. |sj) reveals one broad pulse per period 
with the pulsed fraction of f p = 6.2 ± 1.0%. The MOS 
light curve extracted at / = /o shows a remarkable sim- 
ilarity with that obtained from the pn data (Fig. ||), al- 
though the MOS time resolution makes its pulse some- 
what broader and pulsed fraction smaller than those in 
the pn light curve. We also note that the ROSAT PSPC 
observation of RX J0806. 4-4123 in October 1993 yielded 



Time (s) 



Fig. 3. EPIC light curves with a time resolution of 60 s 
derived from the 0.12-1.2 keV energy band. The back- 
ground with average values of 3.9-10~ 2 cts s _1 , 3.1-10 -3 
cts s- 1 and 4.5-10" 3 cts s" 1 for pn, MOS1 and MOS2, 
respectively, is not subtracted. Start time is 2000, Nov. 8, 
13:45:12 UT. 



228 degrees of freedom, dof). The low value for the absorp- 
tion may be caused by calibration uncertainties mainly 
below 0.5 keV. The radius of the emitting area derived 
from the blackbody fit is R°° ~ 0.6 (d/lOOpc) km. Fitting 
the pn and MOS spectra separately yields blackbody tem- 
peratures of kT°° = 97 ± 1 cV and 91 ± 1 cV for the pn 
and MOS data, respectively (xl = 2.5 for 224 dof). The 
model fluxes in the 0.1-2.4 keV range, 2.82-10 -12 erg cm~ 2 
S"\ 3.08-10~ 12 erg cm" 2 s" 1 and 3.38-10~ 12 erg cm" 2 s" 1 
as derived from the combined fit for the pn, MOS1 and 
MOS2 spectra (respectively), show differences between the 
three instruments larger than the typical statistical errors 
of 1-2% obtained from the fit. For comparison, Haber] 



signal of such a small pulsed fraction (the signal should 
result in a too low peak with Z\ ~ 1, which cannot be 
discriminated from the noise). 

2.3. The X-ray spectrum of RX J0806. 4-4123 

The X-ray spectra were extracted using single- and 
double-pixel events (pattern 0-4) from pn and pattern 0-12 
events from MOS data. Detector response matrix files were 
used as available from the hardware groups in November 
2001 (version 6.1 for pn and v9q20t5r6 for MOS). First, 
an absorbed blackbody model with normalization param- 
eters allowed to be free for the individual instruments was 
simultaneously fitted to the three spectra (left panels in 
Fig. H). From the best fit a blackbody temperature of 
kT°° = 94±1 eVQ and an upper limit for the hydrogen 
column density nn of 10 18 cm~ 2 were found (xl = 2.9 for 



too sc anty statistics (about 900 counts) for searching a ct al. (1998[ ) derived a flux of (2.9±o 3 g)-10~ 12 erg cm 



s- 1 (90% confidence) from the ROSAT PSPC data. 

Nonmagnetic neut ron star hydrogen atmosphere mod- 
els ( ^avlin et al. 1996 ) ^ fit the pn and MOS spectra fairly 
well, resulting in the surface temperature kT°° ~ 26 eV 
and a 110 pc distance to RX J0806. 4-4123 assuming the 
"canonical" neutron star radius i?NS = 10 km and mass 
M NS = 1.4 M e (xl = 1-8 for 228 dof; right panels in 
Fig. ^|). To account for systematic calibration uncertain- 
ties, nn was allowed to vary in fits to the spectra from the 
different instruments, while the surface temperature and 
distance were kept at the same values for all three spec- 
tra. The obtained values of the hydrogen column density 
are n H ~ 1.6-10 20 cur 2 , 1.0-10 20 cm" 2 and 0.4-10 20 cm" 2 
for the fits to the pn, MOS1 and MOS2 data, respectively. 
Neither pure iron nor solar-mixture atmosphere models fit 
the EPIC spectra of RXJ0806. 4-4123 {xl ^ 20). 



2 The superscript °° is used to indicate that quantities are 
given as measured by a distant observer. Errors on the spectral 
parameters are given at a 90% confidence level. 



The hydrogen atmosphere models are ava ilable at 
tittp : //legacy .gsfc .nasa.gov/docs/xanadu/xspec 
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Fig. 5. Pulse profiles extracted from the pn (top 
panel) and combined MOS1 and MOS2 data from 
RXJ0806. 4-4123 at / = / = 0.087940 Hz. 



For instance, if a neutron star is rather massive, so that 
[M N s/M ]/[i? NS /l km] > 0.177, the whole star's surface 



Fig. 4. Power spectrum (top panel) and frequency depen- 
dence of the Odds-ratio. The dashed line in the top panel 
shows a Z\ = 21.5 level below which all peaks can be 
regarded as noise with a probability of 100%. 



3. Discussion 

The XMM-Newton observations of the INS candidate 
RX J0806. 4-4123 revealed pulsations in its X-ray flux with 
a period of 11.3714 s. This is well within the range of 5.16 
s to 22.7 s observed from three other dim INSs. Although 



detecte d at a statistical level of 3.5g, the similar pulse 



becomes visible (e. g. Zavlin et al. 1995b ), which should 
result in a strong smearing of the modulated X-ray flux 
emitted by the neutron star (even if the radiation origi- 
nates from a small spot on the star's surface). However, 
assuming that radii and masses of the detected INSs are 
clustering around the canonical values of i?NS = 10 km 
and Mns = 1.4 Mq, we may conclude that the different 
observed pulsed fractions are mainly caused by geometri- 
cal effects, namely, different orientations of the stars' axes 
and/or different sizes of the emitting areas. 

In this respect it is interesting that for the two dim 
INSs with the estimates on strength of the s urface mag- 
netic fields, -2-10 13 G for RX J0720. 4-3125 flZane et al 
2002 1 ) and -2-10 14 G for RBS 1223 flHambaryan et al 



profiles seen in the data from all EPIC instruments in- 
spire confidence that the modulation is real and indicates 
the rotation period of the neutron star. The pulsed frac- 
tion of 6.2% is the weakest detected so far in X-rays from 
these objects (9% in the X-ray flux of RX J0720.4-3125, 



Cropper et al. (200l| ); 20±2% in RBS 1223 [Hambaryan 
et al (|2002D ; 43±14% in RX J0420. 0-5022, |Haberl et al 
(1999(2 



2002 ), the pulsed fraction in the X-ray flux from the object 
with the stronger field is significantly higher, as is expected 
if the radiation is emitted from hot spots around the mag- 
netic poles. Hot polar caps can either be produced by a 
flux of relativistic particles falling down onto the polar 
caps from the pulsar magnetosphere or anisotropy of the 
thermal conductivity in the neutron star crust caused by a 



Such a modulation is too small to be found in very strong magnetic field (see, e. g., Greenstein & Hartke 



available data of any of the other dim INSs. For example, 



only an upper limit of 4.5% was derived from a 450 ks 



Chandra observation of RX J1856. 5-3754 (Ransom et al 
2002|)~| A simplest explanation for a small pulsed frac- 
tion in X-ray radiation from an INS is that the rota- 
tional axis of the star is nearly co-aligned with either the 
line of sight to a distant observer or the magnetic axis of 
the star (or both). Another reason resulting in the small 
pulsed fraction can be the effect of light bending in the 
strong gravitational field near the neutron star surface 
(e. g. ^echenick et al. 1983j; [Riffert fc Meszaros 1988Q, 



1983). While it looks feasible to explain the shallow mod- 
ulation observed from RXJ0806. 4-4123 by the anisotropic 
heat conduction, Cropper et al. (200l| ) argue that the 
relatively smooth temperature gradient on the neutron 
star surface produced by this mechanism alone can only 
marginally explain the pulse profile of RXJ0720. 4-3125. 
Higher pulsed fractions will be even more difficult to rec- 
oncile with this model, unless the surface radiation is in- 
trinsically anisotropic as it m ay emerge from a magne- 
tized neutron star atmosphere ( Pavlov et al. 1994 ; Zavlin 
et al. 1995a). It is clear that whatever mechanism pro- 
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Fig. 6. EPIC spectra of RXJ0806. 4-4123 together with the best fit blackbody (left) and atmosphere (right) models. 
Most of the residuals (shown in the bottom panels) are probably caused by calibration uncertainties around instru- 
mental absorption edges. 



duces the hot polar regions, measuring rotational period 
changes and, therefore, magnetic field strengths for the 
pulsating dim INS is crucial for understanding the origin 
of their X-ray emission. 

The EPIC spectral data on RX J0806.4-4123 are well 
consistent with a smooth, featureless model and can be 
fitted with both blackbody and neutron star hydrogen 
(or helium) atmosphere spectra (since there are still fairly 
large uncertainties in the calibration of the instrument re- 
sponses at energies below 0.5-0.6 keV, we do not consider 
the blackbody approach as rejected by the data because 
of the rather poor quality of the spectral fit). The derived 
blackbody temperature of ~94 eV is within the range of 
values found for other X-ray dim neutron stars, in partic- 



ular, it is close to that of RXJ1605. 3+3249 jMotch et al 
1999|)~|Extrapolating the best-fit blackbody model to the 
optical wavelength range gives a B magnitude tob — 30.8 
(neglecting extinction) for a possible optical counterpart 



- much fainter flux than the upper limit put by Haberl 
et al. |1998| ). It is worthwhile to mention that the opti- 
cal counterparts found for the INSs RXJ1856. 5-3754 and 
RXJ0720. 4-3125 are significantly brighter than the ex- 
trapolations of the blackbody spectra obtained from the 



X-ray data on these two objects (Pons ct al. 2002; Motch 
& Habjrl 1998). The hydrogen atmosphere model fit pre- 



dicts tob — 26.9, still in agreement with the upper limit 
tob > 24. On the other hand, optical fluxes predicted by 
the light-element atmosphere models for RX J1856. 5-3754 



( Pavlov et al. 1996|; |Pons et al. 2002| ) and RX J0720.4-3125 
( Motch fe Haberl 2002[) arc much brighter than the actu- 
ally measured values. If RX J0806. 4-4123 has properties 
similar to those of these two INSs, then the B magnitude 
of its optical flux is expected to be in the range of 28- 
30. Utilizing the improved X-ray position derived from 



the XMM-Newton data on RXJ0806. 4-4123, more strin- 
gent constraints on the neutron star parameters can be 
obtained from future deep optical observations. 
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